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PREFACE 


Retrospectively  the  path  of  progress  seems  exceedingly  clear.  Wrong  turns  seem 
all  but  too  obvious.  During  the  work  on  this  program  a  number  of  previous  interim 
reports  have  been  written.  It  would  appear  to  have  been  a  wise  decision  to  numbers 
these  reports  as  a  series.  This  was  unfortunately  not  done  originally.  Nevertheless, 
hoping  not  to  seem  presumptuous  the  following  numbering  will  be  used  subsequently: 

Interim  Report-I  —  NASA  Contractor  Report  3670  "A  Study  of  the  Stress  Wave 
Factor  Technique  for  Characterization  of  Composite  Materials",  E.  G.  Henneke,  II,  J. 
C.  Duke,  Jr.,  W,  W,  Stinchcomb,  A.  Govada,  and  A.  Lemascon  (1983). 

Interim  Report-II  —  NASA  Contractor  Report  174870  "A  Study  of  the  Stress 
Wave  Factor  Technique  for  Characterization  of  Composite  Materials",  A.  K.  Govada, 
J.  C.  Duke,  Jr.,  E.  G.  Henneke,  II,  and  W.  W.  Stinchcomb  (1985). 

Interim  Report-Ill  —  NASA  Contractor  Report  4002  "A  Study  of  the  Stress 
Wave  Factor  Technique  for  Characterization  of  Composite  Materials"  A.  Sarrafzadeh 
(-Khoee),  M.  T.  Euernan,  J.  C.  Duke,  Jr.,  and  E.  G.  Henneke,  II  (1986). 

In  addition  to  instituting  the  numbering  sequence  the  introduction  of  this  report 
is  meant  to  serve  as  a  state-of-the-art  review  of  our  understanding  of  the  acousto- 
ultrasonic  method  as  applied  to  the  nondestructive  evaluation  of  advanced  materials. 
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INTRODUCTION 


Background 

The  original  intention  of  Alex  Vary,  the  originator  of  the  acousto-ultrasonic 
(AU)  method,  was  to  provide  a  measurement  procedure  that  would  be  able  to  quan¬ 
tify  the  integrated  effect  of  the  globally  distributed  damage  characteristic  of  fiber 
reinforced  composite  materials. ^  The  original  evaluation  procedure  provided  a 
"stress  wave  factor"  that  correlated  closely  with  several  material  performance  parame¬ 
ters.  The  original  procedure  involved  introducing  a  mechanical  excitation  into  a 
sample  by  means  of  a  piezoelectric  ultrasonic  transducer  and  detecting  the  resulting 
disturbance  at  another  point  a  short  distance  (1  in.,  3  cm)  away  on  the  same  surface 
using  an  acoustic  emission  sensor.  This  received  voltage  signal  was  directed  to  a 
counter  that  counted  the  oscillations  above  a  selected  threshold  voltage  per  selected 
time  interval,  while  the  excitation  occurred  repetitively.  This  quantification  proce¬ 
dure  was  examined  in  detail  by  Henneke  et.  al.  in  NASA  CR  3670  (Interim  Report-I). 

Since  the  introduction  of  this  evaluation  procedure  it  has  been  considered  for 
application  to  a  variety  of  situations,  including: 

-  advanced  composites 

-  bonding  of  hybrid  structures 

-  adhesive  bond  strength 

-  wood  products 

-  wire  rope 

Although  the  original  approach  had  as  its  objective  materials  characterization 
the  procedure  has  been  adapted  for  flaw  detection  in  some  instances.  The  activity  of 
the  research  program  reported  here  has  focused  primarily  on  the  development  of  a 
fundamental  understanding  and  the  advancement  of  fundamental  applications  that 
utilize  the  AU  method  for  characterization  of  materials  for  which  detection,  location 
and  identification  of  imperfection  can  not  at  present  be  analyzed  with  satisfactory 
mechanical  performance  prediction  models.  This  is  not  to  say  that  the  use  of  AU  for 
flaw  detection  is  unimportant,  but  rather  to  point  out  that  the  work  reported  here  has 
not  been  directed  at  optimizing  this  capability;  the  opposite  may  at  time  be  the  case. 

Now  as  with  all  nondestructive  evaluation  method  development  it  is  essential  to 
establish  as  clearly  as  possible  the  nature  of  the  material  attribute  that  is  the  subject  of 
the  evaluation.  Generally  two  situations  need  be  considered:  On  the  one  hand  is  the 
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question  "What  size  must  the  maximally  severe  imperfection  be  before  detection  is 
possible?",  while  on  the  other  we  wonder  "How  severe  must  the  globally  distributed 
damaged  condition  become  before  detection  is  possible?".  If  accept/reject  criteria  are 
available  the  identification  and  engineering  development  of  the  NDE  method  may  be 
carried  out.  This  is  not  the  case  advanced  composite  materials  and  this  has  made  this 
approach  more  challenging  in  that  the  size  of  the  maximally  severe  imperfection,  e.g. 
a  hole,  and  the  severity  of  globally  distributed  damage  are  easily  within  the  detection 
capabilities  of  NDE.  Therefore  the  objective  -  predicting  the  stiffness  and  durability 
of  a  advanced  material  component  -  is  being  pursued  from  at  least  two  directions. 
One  in  which  modelling  is  directed  at  understanding  how  damage  will  degrade  the 
mechanical  performance  while  the  other  seeks  to  relate  the  nondestructive^  evalu¬ 
ated  mechanical  performance  of  the  degraded  material/component,  using  the  AU 
method  or  other  nondestructive  materials  characterization  methods,  to  the  anticipated 
inservice  performance. 

Physical  Basis-  Philosophically  the  AU  approach  endeavors  to  examine  the  mechani¬ 
cal  performance  of  a  material  using  a  mechanical  excitation  that  propagates  in  a 
direction  similar  to  that  of  inservice  loading.  Since  the  mechanical  excitation  is 
dynamic  its  interaction  with  the  material  and  any  damage  present  is  expected  to  be 
similar  to  the  interaction  that  would  occur  with  stress  released  during  additional 
damage  initiation  or  development.  The  magnitude  of  the  mechanical  excitation 
associated  with  the  AU  method  is  globally  quite  small  and  generally  is  assumed  to 
affect  only  linear  elastic  reactions.  The  actual  effect  in  damaged  regions  is  unknown. 
Whether  it  is  reasonable  to  anticipate  a  correlation  between  performance  at  this  level 
of  mechanical  loading  and  that  present  at  failure  of  the  component  is  a  continuing 
matter  of  consideration  in  this  program. 

Wave  Propagation-  With  the  introduction  of  a  short  duration  mechanical  excitation 
into  a  material  consideration  of  the  nature  of  the  propagation  of  the  energy  associated 
with  the  excitation  away  from  the  point  of  excitation  must  certainly  involve  the  discus¬ 
sion  of  wave  propagation.  The  details  of  such  a  discussion  depend  on  the  assumed,  or 
known  description  of  the  excitation,  the  properties  of  the  material/structure,  and  the 
boundary  conditions.  Efforts  to  explain  aspects  of  the  behavior  observed  using  the  AU 
procedure  by  Williams  et.  al.^  and  Kautz^  have  utilized  ray  analysis,  while  workers  at 
VPISU'^  have  considered  plate  modes;  more  regarding  recent  work  is  included  later  in 
this  report. 
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In  1987  the  Acoustical  Society  of  America  formed  a  committee  on  "Structural 
Acoustics  and  Vibrations"  that  has  focus  consideration  on  the  macroscopic  aspects  of 
acoustic  pressure  disturbances  associated  with  structures.  The  authors  see  significant 
similarities  and  have  coined  the  label  "Microstructural  Acoustics"  to  generally 
describe  work  directed  at  explaining  the  nature  of  the  wave  propagation  and  energy 
partitioning  in  composite  materials  where  the  fiber  architecture/matrix  combination 
constitute  the  structure. 


Techniques  for  Implementation 


The  AU  method  is  very  sensitive  to  the  material  condition  and  as  such  may 
appear  to  provide  results  that  are  difficult  to  reproduce  if  care  is  not  taken.  Figure  1 
is  a  photograph  of  a  full  field  moire  interferometric  pattern  associated  with  the  non- 
uniform  part  of  the  displacement  that  occurs  in  a  graphite  epoxy  laminate  that  has 
been  loaded  in  uniform  uniaxial  tension.  The  nonuniformities  must  be  the  result  of 
imperfections  in  the  microstructural  arrangement  of  fibers  and  matrix,  and  laminae. 


Fig.  1.  Photograph  of  a  moire  interferometric  fringe  pattern  that  provides  a  full-field 
map  of  the  in-plane  axial  displacement.  The  pitch  of  the  reference  grating 
has  been  adjusted  to  subtract  out  the  uniform  axial  tensile  elongation  leaving 
only  the  nonuniform  portion.  (Courtesy  of  the  VPISU  Photomechanics  Labo¬ 
ratory;  D.  Post,  Director) 
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Such  nonuniformities  will  influence  AU  measurements  and  can  be  the  cause  of  the 
lack  of  reproducibility  that  may  be  observed.  Therefore  extreme  care  must  be  taken 
when  positioning  transducers  for  AU  measurement.  If  contact  transducers  are  uti¬ 
lized,  in  addition  to  placement,  the  orientation  of  the  transducer  elements  as  well  as 
the  amount  and  distribution  of  the  impedance  matching  couplant  must  be  repeated  to 
achieve  reproducibility.  If  composite  materials  containing  surface  breaking  damage 
are  being  considered  noncontact  or  dry  coupled  transducers  are  recommended.  Wet 
coupled  transducers  may  be  used  successfully  if  the  dwell  time  of  the  couplant  is  short. 
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Fig.  2.  Generalized  schematic  of  the  Acousto-Ultrasonic  (AU)  method 


Figure  2  is  a  generalized  schematic  of  the  AU  measurement  arrangement. 
Depending  on  the  material  system,  the  measurement  objective,  and  the  application 
constraints  numerous  possible  combinations  of  source,  detector(s)  and  data  analysis 
are  possible. 
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Instrumentation 


Transducers  -  Figures  3-6  are  examples  of  several  source  and  detector  combinations; 
both  voltage  versus  time  and  amplitude  spectrum  plots  are  provided.  A  Panametrics 
Model  5052UA  pulser  receiver  and,  where  required,  a  Panametrics  preamplifier  were 
used  along  with  a  General  Research  PCD  AS  acquisition  board  installed  in  an  ARC 
Turbo  XT  to  obtain  the  results.  Figure  3  a  and  b  are  results  from  using  a  5  MHz,  0.5 
in.  diameter,  (Panametrics  V109)  ultrasonic  transducers  as  both  source  and  detector; 
the  damping  setting  was  2000,  the  energy  setting  was  110  and  the  gain  was  46dB. 
Figure  4  a  and  b  are  results  from  using  a  0.25  in.  diameter,  5  Mhz,  (Ultran  KD25-5) 
dry  coupled  transducer  as  a  detector  and  the  same  source  as  in  Fig.  3;  the  gain  was 
increased  to  66dB.  Figure  5  a  and  b  are  results  from  using  a  pulsed  laser  excitation 
(Kigre  MK-480  YAG  Laser  System)  as  the  source  and  a  100  kHz  Panametrics  AE 
sensor  as  the  detector;  the  preamplifier  was  used  to  provide  40dB  gain.  The  pulsed 
excitation  occurred  on  a  1"  x  6"  Gr/Ep  specimen  4  inches  from  the  detector  location. 
Figure  6  a  and  b  are  results  from  using  the  same  source  as  that  used  in  Fig.  3  and  a 
conical  "Proctor"  transducer  as  a  detector. 


Fig.  3.  a)  Voltage  versus  time  and  b)  amplitude  spectrum  for  the  signal  detected 
using  a  0.5  in.  diameter,  5  MHz  transducer  as  a  source  and  another  as  the 
detector. 
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Fig.  4.  a)  Voltage  versus  time  and  b)  amplitude  spectrum  for  the  signal  detected 

using  a  0.5  in.  diameter,  5  MHz  transducer  as  a  source  and  a  0.25  in.  diameter 
5  MHz  dry  coupled  transducer  as  the  detector. 
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Fig.  5.  a)  Voltage  versus  time  and  b)  amplitude  spectrum  for  the  signal  obtained 

using  a  pulsed  laser  source  and  a  0.5  in  diameter,  100  kHz  AE  transducer  as 
the  detector. 
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Fig.  6.  a)  Voltage  versus  time  and  b)  amplitude  spectrum  for  the  signal  obtained 

using  a  0.5  in  diameter,  5  MHz  transducer  as  the  source  and  a  conical  piezoel¬ 
ectric,  "Proctor"  transducer  as  the  detector. 


Previous  work  by  Rebello  and  Duke^  reported  results  obtained  using  a  pencil 
lead  break  source  and  a  capacitance  sensor  as  a  detector.  For  applications  where  the 
component  is  of  conducting  material  electromagnetic  acoustic  transducers  are  worthy 
of  consideration^.  In  addition  it  is  worth  noting  as  well  that  multiple  transducer  AE 
monitoring  arrangements  are  capable  of  being  used  for  AU  measurement  without 
significant  modification.  In  fact  many  systems  utilize  one  or  more  of  the  sensors  as 
sources  for  transducer/system  calibration.  Such  a  procedure  is  reasonable  for  AE 
monitoring  prior  to  the  occurrence  of  damage,  e.g.  before  beginning  of  service  or 
proof  testing,  or  with  components  that  generally  fail  due  to  propagation  of  localized 
self-similar  damage  development,  i.e.  cracks.  With  components  that  sustain  global 
damage  changes  observed  upon  "calibration"  during  service  or  testing  may  be  the 
result  of  changes  in  the  condition  of  the  material  of  the  component.  Therefore  in 
addition  to  being  able  to  detect  the  occurrence  and  propagation  of  damage  through 
conventional  AE  monitoring  the  "calibration"  procedure  can  be  used  to  perform  AU 
measurements  that  would  indicate  the  integrated  effect  of  the  sources  of  the  AE.  One 
might  imagine  using  AE  source/zone  location  to  identify  areas  for  which  AU 
measurement  would  then  be  accomplished  by  using  the  appropriate  AE  transducers  as 
sources  with  others  acting  as  AU  detectors. 


Introduction 


7 


Data  Acquisition/ Analysis-  The  approach  to  acquiring  and  analyzing  the  AU  data  for 
engineering  application  or  research  development  may  differ.  Generally  for  research 
the  emphasis  is  placed  on  flexibility  while  with  the  engineering  application  speed  of 
analysis  is  at  a  premium.  In  the  Materials  Response  Laboratory  at  VPISU  the  over¬ 
riding  objective  is  to  understand  the  mechanical  performance  of  advanced  material 
systems,  so  consequently  it  is  desirable  to  consider  the  observed  AU  response  from  a 
variety  of  perspectives.  However,  once  the  analysis  procedure  is  selected  and  refined, 
by  stripping  flexibility,  considerable  economy  as  regards  time  to  analyze  the  AU 
performance  is  possible. 

The  system  used  for  most  of  the  results  presented  in  this  report  includes  a  PC 
based  digital  data  acquisition  board,  either  the  General  Research  1  channel,  20  MHz 
real-time  (160  MHz  equivalent  time)  sampling,  4096  point  record  length  PCD  AS,  or 
the  Sonotek  2  channel,  25  MHz  real-time  (200  MHz  equivalent  time)  sampling,  64k 
record  length  STR*825,  and  analysis  software  developed  in-house.  It  should  be 
pointed  out  that  analog  data  analysis  is  preferable  in  some  applications,  but  digital 
data  storage  and  analysis  is  significantly  more  convenient  during  the  development 
phase. 

The  analysis  software  developed  has  been  described  in  some  detail  earlier  by 
Kiernan  and  Duke^.  Emphasis  has  been  placed  on  considering  the  spectral  domain 
results  in  an  effort  to  identify  various  modes  of  energy  propagation.  The  component 
geometry,  the  matrix/fiber  architecture  of  the  material  and  the  damage  condition  are 
all  seen  as  contributing  factors  to  how  the  energy  is  partitioned  into  various  modes  of 
propagation.  By  considering  the  nature  of  the  source,  the  boundary  conditions  and 
the  material  properties  progress  is  being  made  in  unfolding  the  complicated  collection 
of  modes  participating.  Identifying  the  modes  is  expected  to  lead  to  a  more  complete 
understanding  of  why  the  AU  behavior  correlates  with  mechanical  performance.  In 
the  meantime  as  this  understanding  increases  the  exploration  of  the  correlation  with 
various  types  of  damage  that  occurs  in  advanced  composites  continues.  To  simplify 
the  chore  of  comparing  the  AU  response  from  point  to  point,  or  from  material  to 
material  the  use  of  a  method  developed  to  quantify  random  noise^  was  suggested.^ 
The  purpose  was  to  quantify  the  nature  of  the  AU  signal  by  means  of  a  small  set  of 
parameters  or  stress  wave  factors  making  comparison  much  easier  than  that  required 
to  compare  the  entire  temporal  or  spectral  domain  information.  This  approach  def¬ 
ines  several  quantities  in  terms  of  spectral  moments  that  have  physical  meaning  when 
applied  to  "stationary"  noise  signals.  Since  noise  signals  are  random  these  measures 
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are  essentially  statistical.  Care  must  be  exercised  however  in  translating  the  original 
interpretation  over  to  the  AU  signals. 


M 


n  “ 


’  f 
■^max 

0 


f"S(f)df 


Mq  is  the  zeroth  moment  of  the  Power  Spectral  Density  Distribution  and  is 
proportional  to  the  energy  of  the  signal.  M;[/Mo  is  the  centroid  frequency  of  the 
power  spectrum. 


The  observation  of  favored  energy  bands  in  the  spectral  domain  of  AU  signals  has 
motivated  a  specialization  of  this  method. 


Al(fa,fb) 


S(f)df 


A2  -  Ml/Mo 


(a  and  b  indicate  upper  and  lower  limits  of  the  frequency  band  of  interest) 


Note:  A’s  are  used  to  denote  versions  of  Rice’s  original  measures  that  have  been 
specialized  for  purposes  of  AU  evaluation. 


By  evaluating  the  area  under  a  portion  of  the  power  spectral  density  distribution 
of  the  AU  signal  it  is  possible  to  estimate  the  energy  associated  with  particular  modes 
or  combinations  of  modes  of  energy  propagation.  Tracking  the  manner  in  which  this 
measure  changes  provides  greater  sensitivity  to  stiffness  and  durability  degrading 
damage,  than  other  AU  measures  considered. 


The  software  provides  other  means  of  comparing  AU  response  (signals)  as  well 
as  allowing  for  statistical  evaluation  of  a  set  of  AU  response  measurements  in  order  to 
determine  a  single  representative  value. 


Data  Interpretation 

Consider  a  thin  graphite  fiber  reinforced  epoxy  laminate  where  a  5  Mhz,  0.5  in. 
diameter  transducer  is  used  as  a  source  and  a  5  MHz,  0.5  in.  diameter  transducer 
placed  one  inch  away  on  the  same  surface  where  the  line  passing  through  the  trans¬ 
ducer  centers  is  parallel  to  a  fiber  direction.  Figure  7  displays  the  variation  of  Al,  for 
such  an  experimental  arrangement,  as  a  function  of  azimuthal  angle.  The  center 
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position  of  the  source  transducer  is  maintained  in  the  same  location.  The  transducer 
source/detection  tandem  is  however  removed  and  replaced  at  the  various  angles  for 
which  measurements  are  made.  Insets  show  the  spectra  for  the  0®,  45®  and  90® 
orientations.  When  considering  the  spectrum  of  the  signals  detected  in  this  manner 
two  frequency  "peaks",  or  bands  appear  to  contain  a  significant  amount  of  the  energy. 
The  frequencies  MHz  and  «=3.2MHz  are  significantly  less  than  the  center 
frequency  of  the  transducers  used.  The  first  peak  is  consistent  with  the  frequency  of 
the  fundamental  extensional  (symmetric)  plate  mode  and  the  second  with  the  flexural 
(antisymmetric)  mode;  see  section  on  "Microstructural  Acoustics"  and  the  Appendix 
for  more  details. 

As  a  simple  rule  of  thumb  the  energy  from  the  mechanical  disturbance  tends  to 
travel  most  readily  in  the  fiber  directions.  Figure  8  shows  the  variation  of  the  energy 
as  represented  by  A1  as  a  function  of  azimuthal  angle  for  a  [0/90]s  laminate.  A1  for 
the  0®  direction  is  4  times  that  for  the  45®  direction  while  that  for  the  90®  direction  is 
3  times.  This  energy  dissipation  is  of  course  dynamic  with  the  energy  flux  in  fibers 
and  matrix  material  and  from  ply  to  ply  being  influenced  by  the  intimacy  of  the 
mechanical  bonding  that  exists  between  these  various  entities. 

The  software  mentioned  earlier^  is  able  to  compare  to  signals  either  from 
different  locations,  different  specimens,  or  from  the  same  specimen  at  different  points 
during  the  course  of  its  life.  This  aspect  has  not  yet  been  explored  in  a  substantive 
way  so  consequently  nomenclature  appropriate  for  the  description  of  this  data  has  not 
been  developed. 
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Fig.  7.  A1  versus  azimuthal  angle  for  a  [0/0/0/0]s  Gr/Ep  laminate.  The  measure¬ 
ments  were  made  by  rotating  a  pair  of  0.5  in.  diameter,  5  MHz  transducers 
about  an  axis  through  the  center  of  the  location  of  the  source  transducer. 
(Actually  both  of  the  transducers  are  completely  removed  and  then  placed 
into  the  appropriate  position  each  time  a  measurement  is  made.) 
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Philosophy  of  Approach 

Understanding  the  nature  of  energy  propagation  away  from  the  source  of  a 
mechanical  disturbance  in  a  fiber  reinforced  composite  material  is  complicated  by 
more  than  just  the  orthotropic  or  transversely  isotropic  material  symmetry.  The 
material  is  anisotropic  and  inhomogeneous  by  design  on  the  scale  of  the  fiber,  lami¬ 
nae,  and  possibly  the  component  level  and  also  by  lack  of  precise  processing.  The 
distribution  and  orientation  of  fibers  is  usually  not  ideal,  and  material  and  processing 
variations  lead  to  uncertainty  in  mechanical  properties  for  analysis  of  AU  behavior. 
To  systematically  approach  the  investigation  of  laminated  composite  materials  a 
decision  as  to  how  to  deal  with  coupled  influence  of  geometry,  lamina  properties,  and 
laminae  interaction  must  be  made. 


Transverse  Resonance  Modes 

Initially  efforts  in  this  research  program  have  been  restricted  to  situations  in 
which  plate  mode  propagation  is  readily  possible  with  sources  suggested  earlier. 
Preliminary  to  actually  considering  the  case  of  the  fiber  reinforced  composite  plate 
the  case  of  a  "monolithic  isotropic  plate"  was  examined. 

The  classical  treatment  of  plate  waves  by  Rayleigh  and  later  by  Lamb  have 
considered  the  motion  of  thin  isotropic,  homogeneous  plates.  With  AU  method 
the  short  wavelengths  associated  with  the  high  frequencies  necessitate  utilizing  a  3-D 
elasticity  analysis.  A  full  development  for  the  case  of  an  isotropic  plate,  as  well  as  the 
preliminary  development  appropriate  for  an  anisotropic  plate  is  presented  in  the 
Appendix.  For  means  of  comparison  and  to  provide  a  measure  of  confidence  in  the 
explanations  being  put  forth  the  special  case  of  the  simple  thickness,  transverse 
resonance,  modes  of  a  6061-T6  aluminum.  Table  1,  and  an  AS-4Gr/Pr  288  [OJg 
graphite  fiber  reinforced  epoxy  plate.  Table  2,  is  presented  here. 
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TABLE  1:  Description  of  6061-T6  Aluminum  Plate 

E  =  10  X  106  psi  G  =  3.75  x  10^  psi  p  =  2.536  x  (Ib-s2)/in4 
Plate  thickness  t  =  0.0625  in. 

C22  =  C33  =  (l-v)E/(l  +  v)(l  -  2v)  =  15  X  106  psi 
€55=  E/[2(l  +  v)]  =  G 

TABLE  2:  Description  of  AS-4Gr/PR  288  [0]g  Graphite/Epoxy  Plate 

Ell  =  18.3  X  106  psi  E22  =  1.3  X  106  psi  G12  =  0.66  x  10^  psi 
^12  =  0.32  1^23  ~  0.15  p  =  1.487  x  10"4  (lb-s^)/in4 
Plate  thickness  t  =  0.088  in. 

C33  =  (1  - 1'  12^E22/Eii)E22A^  =  1.34  x  10^  psi 

V-  [(1  +  1/23X1 23 -2i^l2^E22/Eii)] 

C55  =  Gi2 

The  resonant  frequencies  for  the  simple  thickness  modes  may  be  obtained  by  setting 
f  =  0  in  eqns.  AlO  -  A32.  The  following  resonant,  or  "cutoff  frequencies  result: 

(1)  u  =  m7rVp/2b  w  =  n7rV5/2b 

(2)  symmetric  case  - 

ui  =  0  U2  =  B‘  sin(m7tX2/2b)  m  =  1,3,5, ... 

U2  =  0  ui  =  Ccos(n7rx2/2b)  n  =  2,4,6, ... 
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(3)  antisymmetric  case  - 

U2  =  0  =  D‘sin(n7rX2/2b)  n  =  1,3,5,  ... 

uj  =  0  U2  =  A‘cos(m7tX2/2b)  m  =  2,4,6, ... 

(4)  Vp  =  (C22/p)-^  Vg  =  (C44/p)-5 


TABLE  3:  Frequencies  and  Displacements  for  Transverse  Resonance  Modes 

Mode  Number _ Symmetric _ Antisymmetric 

0  0  0 


f  =  Vp/4b 

f  =  Vs/4b 

1 

U]^  =  0 

ux  =  D‘ 

U2  =  B‘ 

U2  =  0 

11 

f  =  Vp/2b 

2 

ui  =  C 

u]^  =  0 

U2  =  0 

U2  =  -A‘ 

f  =  3vp/4b 

f  =  3vs/4b 

3 

u^  =  0 

uj  =  -D‘ 

U2  =  -B‘ 

U2  =  0 

f  =  Vs/b 

f  =  vp/b 

4 

ui  =  -C 

uj  =  0 

U2  =  0 

U2  =  A‘ 

The  first  symmetric,  the  second  antisymmetric,  and  the  third  symmetric  modes 
haye  non-zero  displacements  in  the  X2  direction.  Howeyer,  due  to  the  Poisson  effect 
the  modes  with  non-zero  uj  displacements  may  also  conceivably  give  rise  to  a  distur¬ 
bance  in  the  X2  direction. 
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Table  4  provides  the  numerical  values  for  the  aluminum  plate  described  by 
Table  1,  while  Table  5  provides  the  numerical  values  for  the  graphite/epoxy  plate 
described  by  Table  2.  The  Goodier-Bishop  wave  mention  in  the  Appendbc  has  a 
frequency,  fQg  =  2.5  MHz 

Figures  9  and  10  are  amplitude  spectra  of  signals  obtained  from  the  aluminum 
and  graphite/epoxy  plates  respectively.  In  both  cases  a  5  Mhz,  0.5  in.  diameter  trans¬ 
ducer  was  used  as  a  source  and  detector.  The  signal  was  sampled  at  a  rate  of  20  MHz 
and  a  fast  Fourier  transform  of  the  first  512  data  points  provided  the  displayed  spec¬ 
tra.  The  frequencies  corresponding  to  the  various  peaks  have  been  labelled.  Within 
the  bounds  of  experimental  accuracy  a  clear  correspondence  may  be  seen  for  both 
cases. 


Fig.  9.  The  amplitude  spectrum  obtained  using  a  0.5  in.  diameter,  5  MHz  transducer 
as  a  source  and  another  similar  transducer  located  a  center-to-center  distance 
apart  of  1  inch  on  a  6061-T6  aluminum  plate  of  thickness  0.0625  in. 
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TABLE  4:  Cutoff  Frequencies  for  Aluminum  Plate 


6061-T6  (0.0625  inches  thick) 


MODE 

SYMMETRIC 

ANTISYMMETRIC 

NUMBER 

FREQUENCY 

(MHz) 

FREQUENCY 

(MHz) 

0 

0.0 

0.0 

1 

1.95 

.97 

2 

1.95 

3.9 

3 

5.9 

2.9 

4 

3.9 

7.8 

6 

9.8 

4.8 
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.00  2.50  5.00  7.50  10.00 


Fig.  10.  The  amplitude  spectrum  obtained  using  a  0.5  in.  diameter,  5  MHz  transducer 
as  a  source  and  another  similar  transducer  located  a  center-to-center  distance 
apart  of  1  inch  on  a  unidirectional  graphite  fiber  reinforced  plate  of  thickness 
0.088  in. 
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TABLE  5:  Cutoff  Frequencies  for  Gr/Ep  Plate 
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AU  EVALUATION  OF  ORGANIC  COMPOSITE  MATERIALS 


In  turning  attention  to  the  fiber  reinforced  laminates  it  was  decided  to  consider 
as  the  fundamental  element,  or  "critical  element",  after  the  terminology  of  the  structu¬ 
ral  analyst,  a  unidirectionally  oriented  fiber  layer.  The  layer  at  times  is  composed  of 
more  than  one  ply  and  the  orientation  of  the  fibers  in  the  fundamental  element  is 
considered  the  0®  direction.  The  behavior  of  the  fundamental  element  alone  and 
under  various  forms  of  constraint  was  deemed  to  provide  a  systematic  means  of  fully 
assessing  the  AU  and  correspondingly  the  mechanical  performance  of  the  composite 
material.  It  should  be  noted  that  the  means  of  constraining  the  fundamental  element 
are  not  arbitrary  but  are  limited  to  those  situations  that  may  be  developed  experimen¬ 
tally. 


In  addition  to  varying  the  nature  of  the  constraint  of  the  fundamental  element 
through  the  use  of  various  stacking  sequences  additional  conditions  that  result  from 
constraint  degradation  due  to  mechanically  induced  damage  were  considered.  These 
included  matrix  cracking  caused  by  cyclic  loading,  damage  due  to  impact,  cracking  and 
impact  damage  combined,  and  simulated  delamination. 


[O/O/O/O] 

The  first  situation  considered  may  for  consistence  be  considered  the  case  where 
the  two  ply  0°  fundamental  element  is  constrained  by  two  0“  plies.  Figure  11  is  a  plot 
of  the  modulus  and  A1  as  a  function  of  azimuthal  angle;  the  behavior  is  very 
similar.  Figure  12  is  a  plot  of  the  first  element  of  the  transformed  reduced  stiffness 
matrix,  and  Al(3.2,3.4)  as  a  function  of  azimuthal  angle;  Al(3.2,3.4)  drops  off 
more  rapidly  between  20  and  30  degrees  than  does  Qjp  Figure  13  is  a  plot  of  vyi  and 
A2,  the  centroid  frequency,  as  a  function  of  azimuthal  angle;  v  \2  drops  off  more 
sharply  between  40  and  80  degrees  that  does  A2.  Figure  14  is  a  plot  of  the  coefficient 
of  mutual  influence  and  A2  as  a  function  of  azimuthal  angle. 

The  AU  response  before  and  after  central  impact  of  a  3  in.  x  6  in.  plate  was 
measured,  Fig.  15.  The  impact  (1.75  ft-lb,  2.37  joules)  was  caused  by  a  steel  sphere 
approximately  1  in.  in  diameter  dropped  onto  the  specimen  and  captured  upon  the 
first  rebound;  the  plate  was  clamped  along  the  narrow  edges,  Fig.  16.  Damage  due  to 
impact  in  this  unidirectional  laminate  was  splitting. 
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Fig.  11.  Modulus,  Ejj.,  and  A1  versus  azimuthal  angle  for  a  thin  unidirectional 
graphite  fiber  reinforced  epoxy  laminated  plate. 
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Fig.  12.  Qj  j  and  Al(3.2,3.4)  versus  azimuthal  angle  for  a  thin  unidirectional 
graphite  fiber  reinforced  epoxy  laminated  plate. 
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POISSON'S  RATIO  A2  (Thousands) 


Fig.  13.  Poisson  Ratio,  v  and  A2  versus  azimuthal  angle  for  a  thin  unidirectional 
graphite  fiber  reinforced  epoxy  laminated  plate. 
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Fig.  14.  Coefficient  of  mutual  influence  (CMI),  »7xy,x’  ^  versus  azimuthal  angle 

for  a  thin  unidirectional  graphite  fiber  reinforced  epoxy  laminated  plate. 
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[O/O/O/O]  A1 


[0/0/0/01  A2 


Fig.  15.  AU  versus  azimuthal  angle  measured  before  and  after  central  impact  at 
location  of  the  source  transducer  for  a  [O/O/O/O]  Gr/Ep  plate,  a)  A1  unda¬ 
maged  alone,  b)  A2  undamaged  alone, 
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[O/O/O/O]  ftl(.8,1.4) 


[O/O/O/O]  Al(3.2,3.6) 


Fig.  15.  (cont.)  c)  Al(.8,1.4)  undamaged,  alone,  d)  Al(3.2,3.6)  undamaged,  alone. 
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[O/O/O/O]  A1 

_ undaHacfed 

e  aftep  iHpact 


Fig.  15.  (cont.)  e)  A1  undamaged  and  after  central  impact,  together;  all  A1  values  are 
normalized  with  respect  to  largest  value  obtained  before  impact.  The  bold 
lines  indicate  the  A1  values  before  impact  and  the  circles  indicate  the  post 


Fig.  16.  Schematic  diagram  of  the  experimental  arrangement  used  to  carry  out  the 
impact  loading  of  the  specimens;  h  «=  10  feet. 
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[45/0/0/45] 

By  placing  single  45°  plies  on  the  outside  of  the  fundamental  element  it  was 
hoped  to  provide  a  situation  of  constraint  intermediate  to  that  of  the  crossply  configu¬ 
ration  considered  next.  Figure  17  is  a  plot  of  A1  versus  azimuthal  angle. 

A1  (IE-4) 


AZIMUTHAL  ANGLE 


Fig.  17.  A1  versus  azimuthal  angle  for  a  [45/0]s  Gr/Ep  laminate. 
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Work  is  presently  in  progress  to  determine  the  AU  response  before  and  after 
impact  for  this  ply  configuration  and  also  before  and  after  cyclic  loading  to  introduce 
matrix  cracking,  and  after  subsequently  impacting  these  cyclically  loaded  specimens. 


[0/90/90/0]  and  [90/0/0/90] 

By  placing  two  90°  plies  next  to  the  fundamental  element  the  most  extreme 
constraint  in  the  transverse  direction  is  obtained  while  the  constraint  in  the  0°  direc¬ 
tion  is  as  small  as  possible  for  any  neighboring  ply  orientation.  Figure  18  displays  the 
response  of  a  [0/90/90/0]  before  and  after  impact,  and  Fig.  19  is  the  corresponding 
result  for  a  [90/0/0/90]  arrangement.  It  should  be  noted  that  the  value  of  A1  in  the 
direction  that  coincides  with  the  outer  ply  orientation  is  25%  higher  than  for  the 
direction  that  sponds  to  the  internal  ply  fiber  orientation. 

[0/90/90/0]  Al 

_ undanasied 

o  after  iHpact 


o 


Fig.  18.  Al  versus  azimuthal  angle  measured  before  any  mechanical  loading,  bold 
lines,  after  impact,  (circles),  of  a  [0/90/90/0]s  Gr/Ep  specimen. 
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[90/0/0/901  A1 

_ undawacfed 

o  after  iHpact 


Fig.  19.  A1  versus  azimuthal  angle  measured  before,  (bold  lines),  any  mechanical 
loading,  after  impact,  (circles),  of  a  [90/0/0/90]s  Gr/Ep  specimen. 


In  order  to  degrade  the  constraint  of  the  fundamental  element  provided  by  the 
90°  plies  the  specimens  were  cyclically  loaded  in  tension  with  R  =  0.1  in  order  to 
introduce  matrix  cracking.  Evidence  of  the  matrix  cracking  was  confirmed  by  means 
of  the  edge  replication  method.  (Penetrant  enhanced  radiography  had  been  consid¬ 
ered  but  was  not  used  after  it  was  determined  that  the  penetrant  affected  subsequent 
AU  measurements.)  Figure  20  is  a  plot  of  A1  as  a  function  of  azimuthal  angle  typical 
of  a  [0/90/90/0]  graphite  fiber  reinforced  epoxy  laminate  before  cyclic  loading,  after 
cyclic  loading,  and  after  impact  following  cyclic  loading.  A  corresponding  set  for  a 
[90/0/0/90]  is  shown  in  Fig.  21.  Figures  22-25  are  penetrant  enhanced  radiographs  for 
[0/90/90/0]  and  [90/0/0/90]  laminates  that  were  impacted  with  and  without  prior  cyclic 
loading;  1:1  and  2.5:1  versions  are  provided.  Attention  is  drawn  to  the  apparent 
correlation  between  the  magnitude  of  the  A1  measured  before  impact  and  the  extent 
of  damage  development  upon  impact.  Work  is  in  progress  to  further  investigate  this 
phenomenon;  [0/90]3s  and  [90/90/90/0]s  laminates  are  being  studied. 
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[0/90/90/0]  A1 

_  undawaged 

»  aftep  cycling 


[0/90/90/01  A1 

_  after  cycling 

«  after  iwpact 


Fig.  20.  a)  A1  versus  azimuthal  angle  measured  before,  (bold  lines),  any  mechanical 
loading  and  after  cyclic  tensile  loading,  (circles),  and  b)  before  impact,  (bold 
lines),  and  after  impact,  (circles),  of  a  cyclicly  loaded  [0/90/90/0]s  Gr/Ep  spe¬ 
cimen. 
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[90/0/0/901  Al 
_  undawaged 
«  aftep  cycling 


0^ 


Fig.  21.  a)  Al  versus  azimuthal  angle  measured  before,  (bold  lines),  any  mechanical 
loading,  and,  after  cyclic  tensile  loading,  (circles);  and  after  impact,  of  the 
cyclicly  loaded  [90/0/0/90]s  Gr/Ep  specimen  no  measureable  change  in  Al 
was  detected. 


Fig.  22.  Penetrant  enhanced  X-ray  radiograph  of  the  specimen  used  to  obtain  the 

data  shown  in  Fig.  18.  a)  1:1  magnification,  b)  enlargement  of  the  region  near 
the  impact,  2.5:1  magnification. 
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Fig.  23.  Penetrant  enhanced  X-ray  radiograph  of  the  specimen  used  to  obtain  the 

data  shown  in  Fig.  19.  a)  1:1  magnification,  b)  enlargement  of  the  region  near 
the  impact,  2.5:1  magnification. 


Fig.  24.  Penetrant  enhanced  X-ray  radiograph  of  the  specimen  used  to  obtain  the 

data  shown  in  Fig.  20.  a)  1:1  magnification,  b)  enlargement  of  the  region  near 
the  impact,  2.5:1  magnification. 
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Fig.  25.  Penetrant  enhanced  X-ray  radiograph  of  the  specimen  used  to  obtain  the 

data  shown  in  Fig.  21.  a)  1:1  magnification,  b)  enlargement  of  the  region  near 
the  impact,  2.5:1  magnification. 

Because  dynamic  load  redistribution  occurs  during  damage  development  it  is 
believed  that  insight  into  the  nature  of  this  damage  development  is  possible  by  consid¬ 
ering  the  AU  response  of  the  material/component  when  the  stress  waves  caused  by 
the  AU  source  are  somewhat  analogous  to  the  damage  induced  stress  waves. 
Work  is  continuing  on  interpreting  the  observed  AU  response  in  light  of  the  various 
damage  mechanisms  that  may  occur  as  a  result  of  a  particular  mechanical  load  exper¬ 
ience. 
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AU  EVALUATION  OF  CERAMIC  COMPOSITE  MATERIALS 


The  need  for  material  systems  that  will  perform  at  elevated  temperatures  has 
turned  additional  attention  to  the  ceramic  fiber  and  ceramic  matrix  composite  mate¬ 
rial  candidates.  Traditionally  the  toughness  of  monolithic  ceramics  has  been  low 
causing  concern  regarding  their  use  in  designs  that  demand  durable  performance. 
The  inclusion  of  reinforcing  fibers  appears  to  offer  a  means  of  significantly  improving 
the  durability  and  thus  allowing  for  exploitation  of  the  high  temperature  capabilities 
of  these  materials.  Since  with  other  composite  material  systems  the  toughness  attrib¬ 
ute  is  also  associated  with  distributed  damage  development  it  has  been  anticipated 
that  AU  for  characterization  of  ceramic  composite  materials  may  be  advantageous. 

Although  the  shortcomings  of  monolithic  ceramics  are  many,  much  work  is 
being  done  to  develop  "tough"  ceramics.  Toughening  can  be  achieved  in  several  ways, 
including:  the  addition  of  a  second  phase  of  the  same  material,  the  addition  of  par¬ 
ticles  of  another  material  to  form  a  particulate  composite,  or  the  addition  of  whiskers 
or  long  fibers  to  create  a  fibrous  composite  material.  The  roles  of  these  secondary 
materials  are  many,  including: 

1)  To  deflect  an  approaching  crack,  increasing  the  distance  which  it  must  travel 
through  a  specimen,  thereby  increasing  the  life  of  the  material 

2)  To  branch  the  crack,  thereby  eliminating  the  effects  of  one  large  crack  propa¬ 
gating  through  the  material 

3)  To  reduce  the  stress  concentration  at  the  tip  of  the  crack,  thereby  blunting  it 
and  slowing  its  growth 

4)  To  dissipate  fracture  energy  through  fibers  pulling  out  of  the  surrounding 
matrix  material. 

5)  To  prestress  the  material,  with  residual  stress  remaining  after  processing. 

These  composite  materials  have  shown  a  much  greater  strain  energy  to  failure 
compared  to  their  monolithic  counterparts.  In  many  cases  material  properties,  such  as 
strain  to  failure  and  ultimate  strength,  have  been  increased  significantly. 


Literature  Survey 

Imperfection  Detection-  Unfortunately  for  the  researcher,  very  little  information  has 
appeared  in  the  literature  regarding  NDT  of  ceramic  matrix  composite  materials.  It 
seems  to  mean  that  either  research  in  this  area  is  badly  lacking  or  the  work  that  has 
been  done  is  considered  proprietary.  Most  of  the  information  must,  therefore,  be 
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drawn  from  the  studies  of  monolithic  ceramics  and  then  must  be  applied  to  ceramic 
matrix  composite  materials. 

The  available  literature  leads  to  the  area  of  flaw  detection  and  critical  flaw  sizes 
for  a  monolithic  ceramic  materials.  The  size  of  the  largest  imperfection  which  can 
exist  in  a  component  and  not  cause  catastrophic  failure  is  considered  to  be  the  critical 
flaw  size.  Much  of  the  NDT  research  has,  therefore,  been  centered  on  locating 
imperfections  and  determining  the  smallest  flaws  which  can  be  accurately  found  by 
various  NDT  techniques. 

The  purpose  of  using  reinforcement  in  ceramic  materials  is,  of  course,  to  reduce 
sensitivity  of  the  material  performance  to  small  imperfections.  One  goal  of  the 
researcher  should  be  to  identify  the  smallest  strength  controlling  flaw  which  may  exist 
in  a  particular  material.  Detection  techniques  should  then  be  developed  to  locate  any 
imperfections  larger  than  this  size. 

The  critical  flaw  size  in  some  monolithic  materials  has  been  reported  to  be  in 
the  10  micrometer  range.  NDT  techniques  applied  to  monolithic  ceramics  have 
been  reported  which  are  able  to  locate  flaws  in  the  tens  of  micrometers  range,  so 
these  techniques  should  be  easily  applied  to  ceramic  composites. 

One  problem  of  concern  is  that  most  studies  in  this  area  use  imperfections 
imbedded  in  samples  by  the  investigators.  The  usefulness  of  the  techniques  becomes 
questionable  when  applied  to  real  ceramic  components,  out  of  laboratory  conditions. 
Another  problem  is  the  time  required  to  perform  the  inspection;  it  may  take  hours  to 
fully  inspect  even  a  small  region  for  microscopic  imperfections. 

Techniques  which  are  of  use  in  locating  and  sizing  of  small  imperfections 
include  high  frequency  ultrasonic  scans,  ultrasonic  scattering  and  attenuation 
measurements,  and  microfocus  radiography. 

Porosity  and  Density-  Other  imperfections  which  are  desirable  to  locate  in  ceramic 
materials  are  large  pores  and  density  variations,  that  would  have  a  negative  effect  on 
the  strength  or  toughness  of  the  material. 

Most  pores  and  density  variations  result  from  numerous  process  variables.  New 
processes  for  the  production  of  ceramic  matrix  composites  are  being  developed  which 
lower  density  variations  and  porosity. 
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One  technique  being  developed  for  the  production  of  ceramic  matrix  composi¬ 
tes  is  Chemical  Vapor  Infiltration  (CVI).l^  In  this  technique,  a  chemical  vapor  con¬ 
taining  molecules  of  the  matrix  material  is  passed  through  a  preformed  mesh  of  fibers. 
In  the  presence  of  heat  and  pressure,  the  matrix  is  deposited  uniformly  on  the  ceramic 
fibers.  The  process  takes  hours  or  days  to  complete,  but  the  resulting  matrix  is  signifi¬ 
cantly  more  uniform  and  dense  than  composites  produced  by  other  methods. 

Composites  produced  by  CVI  have  been  reported  to  have  theoretical  densities 
as  high  as  85-90%,20»21  however,  a  large  number  of  pores  may  be  distributed 
throughout  the  composite.  Oak  Ridge  National  Laboratories  reports  that  pore  sizes 
as  large  as  several  hundred  micrometers  exist  within  these  composites.  It  has  been 
hypothesized  that  further  densification  will  have  little  beneficial  effect  on  strength 
and  toughness  of  these  materials.^^ 

Several  other  processing  methods  are  mentioned  in  the  available  literature 
including  reaction  bonding,  polymer  precursor  synthesis,  and  sol-gel  techniques,^^ 
however,  the  prevalent  processing  method  is  hot  pressing.  It  is  difficult  to  produce 
uniform  material  using  this  method;  large  variations  in  the  material  density  after 
processing  have  been  reported.  These  variation  will  have  a  pronounced  effect  on  the 
properties  of  the  finished  material.  Such  inhomogeneities  may  be  produced  by  any 
number  of  processing  problems  (such  as  irregularities  in  raw  materials,  hot  pressing 
variations  such  as  time,  temperature  and  pressure). 

Large  pores  often  remain  in  the  material  even  after  hot  pressing.  One  reason  for 
this  is  because  the  reinforcement  fibers  inhibit  the  flow  of  viscous  glass  (or  other 
matrix  material).  The  flowing  matrix  material  may  be  restrained  from  flowing  into 
certain  areas  of  the  mesh,  thereby  leaving  voids  and  pores. 

NDT  techniques  have  proven  to  be  useful  in  several  ways  in  this  area.  Location 
of  pores  larger  than  a  critical  flaw  size  would  be  necessary  to  assure  that  a  component 
would  be  trustworthy  in  service.  High  frequency  ultrasonic  techniques  and  microfocus 
radiography  could  prove  to  be  effect  in  locating  and  sizing  such  flaws. 

Location  and  characterization  of  clusters  of  small  pores  or  density  variations 
would  probably  be  more  useful  in  the  investigation  of  ceramic  composite  materials 
because  of  the  decreased  dependency  upon  single  imperfections.  Ultrasonic 
ing  and  attenuation  measurements  have  been  used  to  measure  density  gradients  in  the 
monolithic  ceramic  materials.  Ultrasonic  velocity  measurements  have  been  applied  to 
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monolithic  ceramics  and  material  densities  variations  as  small  as  2%  have  been  identi¬ 
fied.  X-ray  techniques  have  also  been  used.  These  techniques  may  also  be  applicable 
to  ceramic  matrix  composite  materials.^^’^^’^^’^^ 

Matrix  Variation-  Most  matrix  materials  for  ceramic  matrix  composites  have  been 
ceramics  which  have  been  produced  in  monolithic  form  at  one  time,  however, 
required  toughening  or  strengthening.  The  condition  of  the  matrix  material  is  of  as 
great  an  importance  in  composite  materials  as  the  condition  of  the  reinforcing  fibers. 

Variations  in  grain  size  may  have  an  undesirable  effect  on  the  properties  of  a 
ceramic.  An  oversized  grain  may  be  the  source  of  a  crack  or  if  too  large,  the  grain 
may  exceed  the  critical  flaw  size. 

The  existence  of  different  phases  of  the  matrix  ceramic  can  greatly  change  the 
properties  of  the  resultant  composite.  The  existence  of  secondary  phases  of  the 
ceramic  matrix  could  have  beneficial  or  undesirable  effects  on  the  material.  In  some 
materials  inclusions  of  a  secondary  phase  are  considered  to  be  the  toughening  agent  - 
the  inclusions  deflect  approaching  cracks.  The  secondary  phase  may  also  degrade  the 
overall  performance.  Secondary  phase  inclusions  may  be  crack  initiation  points; 
microcracks  may  occur  in  the  two  phase  matrix  at  a  lower  load  than  would  be  expected 
in  single  phase  matrices. 

Residual  Stress-  As  in  other  composite  systems,  residual  stresses  are  present  in  cer¬ 
amic  matrix  composite  materials  due  to  the  unequal  coefficients  of  thermal  expansion 
between  the  various  components  of  the  composite  system.  The  materials  are  pro¬ 
cessed  at  a  high  temperature,  and  after  cooling,  the  constituents  of  the  composite 
system  attempt  to  shrink  different  amount,  therefore,  creating  residual  stresses. 

Residual  stresses  may  also  be  present  due  to  other  processing  techniques,  such 
as  rolling  or  pressing  of  the  material.  Machining  of  components  may  also  leave 
stresses  in  the  material.  If  the  residual  stress  levels  within  the  composite  are  large 
enough,  microcracking  of  the  matrix  will  result. 

An  understanding  of  these  residual  stress  levels  will  be  important  in  the  design 
and  analysis  of  structures  constructed  of  advanced  composite  materials.  A  measure¬ 
ment  of  these  stress  levels  would  prove  to  be  of  great  use  when  determining  the 
suitability  of  a  part  for  service. 
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Anisotropy-  Both  whisker  and  fiber  reinforced  (long  fibers)  composite  materials  will 
be  anisotropic  to  a  certain  degree.  On  the  other  hand,  randomly  oriented  particulate 
composites  can  usually  be  treated  theoretically  as  statistically  isotropic.  The  variation 
of  mechanical  properties  in  various  directions  may  result  from  fiber  orientation,  hot 
pressing,  or  rolling  of  the  ceramic  composite  during  preparation. 

An  understanding  of  the  degree  of  anisotropy  is  required  for  proper  design  and 
use  of  a  material  system.  A  measure  of  the  relative  material  properties  in  various 
directions  would  be  of  great  use.  This  measurement  may  prove  to  be  difficult  to 
accurately  determine.  A  NDT  method  to  determine  these  properties,  or  at  least  give 
a  qualitative  assessment  of  their  variation,  would  be  useful. 

Presently,  mechanical  testing  is  the  dominant  method  of  determining  material 
properties  in  anisotropic  specimens.  Ultrasonic  velocity  measurement  have  been  used 
on  ceramic  materials  to  determine  anisotropy.^^’^^  Acousto-Ultrasonic  measure¬ 
ments  have  been  shown  to  have  a  direct  correlation  to  direction  in  other  anisotropic 
materials.^® 

Fiber  Distribution-  Homogeneous  fiber  distribution  throughout  a  matrix  material  is 
an  important  quality  in  a  composite  material  system.  It  is  not  uncommon  to  have 
areas  of  fiber  or  whisker  concentration  or  sparsity  in  these  materials.  These  regions 
would  obviously  have  detrimental  effects  on  the  performance  of  the  material.  The 
location  and  measurement  of  degree  of  concentration  variation  could  be  used  in 
quality  control  systems.  A  measurement  of  the  effect  of  these  inhomogeneities  on  the 
overall  material  characteristics  would  be  of  even  greater  use. 

Ultrasonic  velocity  tests  and  attenuation  measurements  could  be  useful  in  locat¬ 
ing  these  imperfections.  X-ray  methods  would  also  be  useful  in  locating  and  sizing  the 
regions  of  concern.  Yet,  with  all  of  these  methods  it  left  to  the  NDT  inspector  to 
determine  the  effect  these  regions  may  have  on  the  performance  of  the  material.  AU 
would  yield  a  interactive  measurement  of  the  material  characteristics  and  could  be 
used  to  determine  if  the  material  was  suitable  for  use. 

Fiber  Degradation-  Several  fibers  used  for  reinforcement  of  a  composite  systems  have 
been  shown  to  degrade  in  the  presence  of  high  temperatures  and  oxidizing 
environments.  These  factors  may  exist  during  material  processing  or  while  in 
service. 
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The  matrix  material  of  a  ceramic  composite  material  generally  protects  the 
fibers  from  the  surrounding  atmosphere.  However,  in  service,  microcracks  may  allow 
exposure  to  degrading  atmospheres. 

For  example,  it  has  been  shown  that  SiC  fibers  may  begin  degrading  at  tempera¬ 
tures  as  low  as  1000°  C.  An  oxide  scale  forms  on  the  fiber  surface  and  the  fiber 
becomes  extremely  notch  sensitive.  After  such  exposure,  fiber  strength  may  be 
reduced  by  as  much  as  75%.^^  The  degradation  of  the  fibers  in  association  with 
matrix  microcracking  will  greatly  reduce  the  stiffness  and  strength  of  a  material. 

It  is  questionable  if  X-ray  methods  would  be  able  to  show  fiber  degradation, 
however,  fiber  breakage  could  probably  be  noticed.  A  reduction  in  stiffness  and 
strength  would  be  found  if  the  part  were  loaded  or  proof  tested,  however,  these 
methods  may  advance  damage  which  already  exists  within  the  structure.  AU  could 
possibly  show  the  extend  of  fiber  degradation  as  well  as  the  resulting  material  charac¬ 
teristics. 

Damage  Development-  A  ceramic  matrix  which  has  been  reinforced  with  continuous 
fibers  not  only  has  decreased  sensitivity  to  microscopic  flaws,  but  also  has  the  ability 
to  sustain  damage  without  catastrophic  structural  failure.  It  has  been  shown  that 
ceramic  matrix  composite  materials  with  long  fibers  for  reinforcement  have  a  general 
failure  pattern:  elastic  behavior  until  the  formation  of  microcracks  within  the  matrix, 
increasing  strain  levels  while  microcracking  appears  throughout  the  matrix,  and  fail¬ 
ure  with  a  single  predominant  crack  through  the  thickness  of  the  specimen.^^’^^’^^ 

The  detection  of  the  onset  of  non-elastic  behavior  would  be  of  great  use  to  the 
designer.  Acoustic  emission  monitored  during  testing  may  be  used  to  locate  the  load 
at  which  the  specimen  begins  to  crack.  This  information  would  allow  for  parts  to  be 
design  to  stay  below  the  load  levels  that  produce  such  cracks. 

This  is  of  little  use  when  considering  composite  components  which  have  been 
used  in  service  or  those  which  may  have  been  exposed  to  fatigue  loading.  Fatigue  type 
loads  would  tend  to  produce  microcracks  within  the  matrix,  however,  a  further  study 
of  the  fatigue  damage  development  is  needed. 

A  measure  of  the  material  integrity  after  some  load  history  is  necessary.  Such  a 
measure  may  be  made  nondestructively  on  ceramic  matrix  composite  materials 
through  the  application  of  AU.  The  stress  wave  which  travels  through  the  composite 
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specimen  interacts  with  the  material  and  any  damage  which  may  be  present,  thereby 
resulting  in  a  measure  of  the  material’s  performance.  This  quantitative  measure  of  the 
material  characteristics  would  eliminate  the  dependence  upon  subjective  decision 
required  when  other  NDT  measurement  techniques  are  use. 

Fiber-Matrix  Interface-  It  has  been  stated  that  the  behavior  of  ceramic  matrix  compo¬ 
site  material  is  highly  dependent  upon  the  characteristics  of  the  fiber  matrix  inter- 

In  a  composite  in  which  the  constituents  are  of  the  same  com¬ 
position,  any  cracks  which  propagate  through  the  material  tend  to  "see"  only  one 
material,  and  therefore,  the  material  acts  as  if  it  were  monolithic.  These  materials 
generally  exhibit  very  brittle  behavior  and  show  little  tolerance  to  damage.  A  compo¬ 
site  material  in  which  a  strong  "bond"  between  the  fibers  and  the  matrix  exist  will  also 
result  in  a  very  brittle  composite  -  the  fibers  will  break,  rather  than  pull  out  of  the 
matrix.  These  composites  offer  little  or  no  resistance  to  a  propagating  crack.  A  fric¬ 
tional  interface  results  in  a  better  performing  composite.  The  fibers  can  pull  out  or 
deflect  approaching  cracks,  thereby  increasing  the  energy  needed  to  fracture  the 
material. 

Several  methods  of  altering  the  fiber  matrix  interface  to  produce  a  better  per¬ 
forming  material  are  used.  The  most  successful  method  seems  to  be  fiber  coatings. 
The  fibers  are  coated  with  a  thin  layer  of  another  material  which  creates  a  frictional 
bond  with  the  matrix.  Because  the  fiber-matrix  interface  has  a  pronounced  effect  on 
the  performance  of  the  composite,  this  area  deserves  much  more  research  attention. 

In-Process  NDT-  It  has  been  reported  that  NDT  has  been  used  to  determine  modifi¬ 
cations  necessary  to  improve  the  quality  and  reliability  of  ceramics  and  ceramic  matrix 
composites.  Various  techniques  have  been  used  including  radiography  and 

ultrasonic  scans. 

Testing  of  specimens  has  been  done  at  various  stages  during  processing,  includ¬ 
ing  NDT  of  raw  materials  (powders  and  fibers),  green  state  ceramics,  and  finished 
materials.  Changes  in  production  that  have  been  suggested  through  nondestructive 
evaluation  have  reduced  density  variations,  porosity,  fiber  misalignment,  and  fiber 
concentrations. 
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Correlation  With  Ultrasonic  C-scan  Results 


Some  fiber  reinforced  borosilicate  glass  matrix  material  has  been  provide,  for 
study,  by  Martin  Marietta.  Prior  to  shipping  these  panels  the  NDE  Group  at  the 
Research  and  Development  Center  of  Martin  Marietta  Corporation  examined  the 
materials  using  ultrasonic  C-scan.  AU  evaluation  of  these  panels  at  VPISU  has 
provided  some  preliminary  results. 

The  panels  were  examined  using  0.5  in.  diameter  5  MHz  transducers  described 
earlier  as  source  and  detector.  The  transducer  pair  was  kept  one  inch  apart  while 
begin  at  one  edge  of  the  4  in.  x  4  in.  plates.  Measurements  were  then  made  at  posi¬ 
tions  every  0.5  in.  travelling  in  a  direction  parallel  to  the  line  through  the  transducers 
centers.  When  a  transducer  reached  the  edge,  the  assembly  was  indexed  perpendicu- 
perpendicular  to  this  line  of  motion,  0.5  in.  and  the  measurements  were  again  made  at 
0.5  in  intervals.  The  signal  was  captured  and  analyzed  in  order  to  determine  the  value 
of  A1  for  each  signal  detected.  The  A1  values  were  assigned  to  the  locations  midway 
between  the  transducers  for  each  location.  Since  the  transducers  were  not  allowed  to 
extend  off  the  plate  and  because  of  their  0.5  in.  diameter  the  AU  measurement  area 
was  3.5  in.  x  2.5  in. 

Figure  26  a  is  a  "topographical"  plot  of  the  measure  A1  for  a  Coming  7740 
borosilicate  glass  plate  reinforced  with  hand  chopped  Nicalon  SiC  fiber  cloth.  Figure 
26  b  is  the  ultrasonic  C-scan  obtained  using  2.25  MHz  transducers  and  squirters. 
Figure  27  a  is  a  "topographical"  plot  of  the  measure  A 1  for  a  Corning  7740  borosilicate 
glass  plate  reinforced  with  hand  chopped  graphite  Fiberite  W-341  cloth.  Figure  27  b 
is  the  ultrasonic  C-scan  obtained  using  2.25  MHz  transducers  and  squirters.  Figure  28 
is  a  wire  frame  plot,  of  the  A1  values  shown  in  Fig.  27.  (A  Corning  7740  borosilicate 
glass  plate  with  6  plies  of  graphite  fiber,  Fiberite  W-341  cloth,  was  also  provided. 
Martin  Marietta  was  unable  to  ultrasonically  scan  the  plate,  and  the  initial  efforts  to 
evaluate  the  AU  response  were  inconclusive.) 

Work  is  continuing  in  an  effort  to  examine  the  correlation  of  AU  behavior  and 
mechanical  performance  for  this  class  of  materials. 
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Fig.  26.  (cont.)  b)  "gray  scale  copy"  of  the  color  normal  incidence  ultrasonic  C-scan 
using  2.25  MHz  transducers  and  squirters.  (Provided  by  Martin  Marietta) 
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Fig.  27.  (cont.)  b)  "gray  scale  copy"  of  the  color  normal  incidence  ultrasonic  C-scan 
using  2.25  MHz  transducers  and  squirters.  (Provided  by  Martin  Marietta) 
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Fig.  28.  "Wire  frame"  plot  presentation  of  A1  values  presented  in  Fig.  27a. 
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LASER  IN/OUT 


Laser  Interferometric  Detection 

In  order  to  apply  the  AU  NDE  method  in  some  situations  noncontact  source 
generation  and  ultrasonic  detection  are  required  because  of  constraints  with  regard  to 
use  of  acoustic  couplant  or  geometric  details  of  the  component.  Although  several 
possibilities  are  available  including  EMATs  and  transducers  coupled  by  means  of  air 
or  elastomeric  materials,  special  consideration  has  been  given  to  pulsed  laser  gener¬ 
ation  of  mechanical  disturbances  and  optical  detection  of  transient  surface  displace¬ 
ments.  Previously  optical  systems  capable  of  measuring  in-plane^^  and  out-of- 
plane^^  displacement  were  developed.  Recently  an  optical  arrangement  that  inte¬ 
grates  these  measurement  procedures  in  order  to  allow  for  making  both  measure¬ 
ments  without  moving  the  specimen  has  been  set  up  in  the  Acousto-Ultrasonics 
Facility  at  VPISU.  Figure  29  is  a  schematic  of  this  optical  arrangement. 


Laser  In 

The  acquisition  of  a  pulsed  laser  system  has  made  it  possible  to  consider  pulsed 
laser  generated  mechanical  excitation  as  an  AU  source.  Figure  5  is  an  example  of  the 
signal  obtained  using  such  a  source  as  detected  by  an  acoustic  emission  sensor.  The 
use  of  such  a  means  of  excitation  however  raises  a  number  of  fundamental  questions 
including: 

1)  What  is  the  nature  of  the  interaction  between  the  laser  generated  excitation 
and  the  material,  both  surface  and  subsurface? 

2)  What  level  of  excitation  is  required  to  provide  a  detectable  signal? 

3)  How  may  the  nature  of  the  excitation  be  controlled  or  modified? 

Work  directed  at  answering  these  questions  is  in  progress. 
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Fig.  29.  Schematic  of  combined  in-plane  and  out-of-plane  transient  displacement 
optical  measurement  arrangement.  The  solid  lines  represent  the  positions 
for  optical  elements  of  the  out-of-plane  detection  setup,  while  the  dashed 
lines  indicate  the  optical  components  used  to  detect  in-plane  displacement. 
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APPENDIX 


The  basic  plate  wave  problem  may  be  formulated  by  considering  the  geometry 
diagrammed. 


The  first  analyses  were  for  long  wavelength  disturbances  in  thin,  isotropic, 
homogeneous  plates: 

Flexural  (Germain-Lagrange) 

A1  Dv2v2u2  +  2/9bU2  =  0  where  D  =  4/ib^/3(l-v) 

Extensional  (Poisson) 

A2  (X'  +  +  U3^3),2  + 

A3  '  +  /i)(U2^2  +  ^3,3)0  /^v2u3  =  /7U3  X 2X/i/(X+^) 


For  these  equations  \x[  -  uj(xi,X3,t). 

For  short  wavelengths,  typical  of  the  AU  method  the  thin  plate  theory  is 
insufficient  since  the  uj’s  depend  on  the  thickness,  therefore  variations  across  the 
thickness  can  not  be  ignored,  ie.  Uj  =  Ui(x2,X2,X3,t). 

A4  +  ^ui,kk  =  P^i 

where 

A5  +  Uj^£) 
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A4  can  be  separated  into  two  uncoupled  systems  of  equations  using: 


A6 

Ui  + 

=  0 

A7 

Vp2v2<^  =  <t> 

CO 

II 

(\+2n)/p 

A8 

Vs2v2Hi  =  Hi 

Vs2  = 

The  conditions  imposed  by  the  boundaries  may  be  described  as: 

A9  02j (xi,±b,X3,t)  =  0 


One  possible  solution  for  the  case  of  SH  waves  finds  ui  =  U2=  0  and 

U3  =  Eiexpi(exi  +  -  wt)  +  E2expi(exi  -  py.2  -  wt)  . 

Another  that  is  of  more  interest  occurs  if  =  H2  =  0  in  eqn.  A6. 

AlO  <p=  <#)q  expi(2CXj[  ±  ax2  -  wt)  =  Vp2  (a2  +  ^2) 

All  H3  =  Hoexpi(^3^  ±  ^X2  -  wt)  w2  =  Vs2(^2  +  ^2) 

A12  6  =  ±tan~^(^/a)  P  wave 

A13  ^2  =  ±tan“^(^//9)  SV  wave 

A  harmonic  wave  motion  in  the  -  direction  with  phase  velocity  v  =  and 
wavelength  X  =  Ztt/C  may  be  seen  to  be  a  superposition,  with  the  appropriate  choice 
of  phase  angles  and  amplitudes,  of  four  waves  of  the  type  described  in  eqn  AlO  and 
four  waves  of  the  type  described  in  eqn  All;  specifically: 

A14  <!>=  (A  sinQ;X2  +  B  cosax2 )  expi  (^x^  -  wt) 

A15  H3  =  i(Csin/5x2  +  Dcos/9x2 )  expi  (^x^  -  cot) 

The  general  form  of  the  displacements  and  stresses  associated  with  these  waves 
could  be  written: 

A16  =  Uj^expi(^X3  “ 

A17  a^j  =  Tijexpi(£xi  -  cot) 
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Now  this  type  of  wave  motion  may  be  divided  into  two  distinct  types: 

1)  symmetric  motion  about  the  xi  -  X3  plane 

A18  Ui  =  i(BCcosaX2  +  003^X2) 

A19  U2  =  (-BasinaX2  +  C^sin/9X2) 

A20  T22  =  -  i?^)cosaX2  +  2CC/9cos^X2] 

A21  T21  =  [-2B^asinQ:X2  +  -  /S2)sin^X2] 

A22  Til  M[B(2q:2  _  ^2  -^2)cosaX2  -  2C^/9cos^X2] 

A23  Fi(a,yS,^)  =  )  2(-;osabsin/9b  +  4a/9^^sinabcos^b  =  0 

2)  antisymmetric  motion  about  the  -  X3  plane 

A24  Ui  =  i(A^sinaX2  "  D/9sin/SX2) 

A25  U2  =  (AQ:COSOiX2  +  0^003/9X2) 

A26  T22  =  “  /92)sinax2  -  2D^^sin^X2] 

A27  T21  =  iju  [  2Aa^  cosq!X2  +  203^X2] 

A28  Til  =  ^  [A(2a2  -  ^2  _  ^2)sinax2  +  2DC^sin^X2] 

A29  F2(o:,/S,£)  ®  "  /32 )  2sina!bcos/9b  +  4a/9C ^sin/ffbco3ab  =  0 

The  two  frequencies  eqns,  can  be  written  in  reduced  form  and  combined  phase 
velocity  calculations  to  yield  the  following  dispersion  relations: 

A30  F(a,/S,C)  =  tan;Sb/tanab  + 

L 

+1  3yininetric,  -1  anti3y]nmetric 
A31  a2  =  (w2/vp2  -  ^2) 

A32  y92  =  (y2/v^2  _  ^2) 
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Extension  to  Anisotropic  Material 

For  the  general  anisotropic  case,  the  following  homogeneous  system  of  equations 
can  be  derived  by  equilibrium  considerations: 


A33  Aj^(C '  '  1  “  ^  ik)  “  0 

where 

C'ijkl  =  Qjkl/^1212  dimensionless  stiffness  matrix  terms 
^  'i  =  2b^i/7r  dimensionless  wave  number 
n  =  w/wg  dimensionless  frequency 
and 

ug  =  ('7t/2b)(Ci2l2/^)^  resonant  frequency  of  shear  wave, 

A]j  is  the  amplitude  of  the  particle  displacement. 

For  a  given  direction  C  f/ 1  ^  1  ,  the  above  expression  yields  three  distinct  roots,  cor¬ 
responding  to  three  different  wave  speeds  and  particle  displacement  directions.  It  is 
also  necessary  to  include  stress  free  conditions  at  the  upper  and  lower  surfaces.  Since 
the  three  waves  are  coupled,  the  situation  is  more  complicated  than  the  isotropic  case. 
The  anisotropic  constitutive  expression  is  significantly  more  cumbersome  making  the 
transcendental  equations  more  awkward  to  handle. 

The  approach  being  followed  in  pursuing  this  analysis  is  as  follows: 

1  Define  a  direction  of  propagation  that  is  of  interest 

2  Formulate  the  equations  of  equilibrium 

3  Describe  the  displacements,  assuming  standing  waves  in  the  X2-direction  with 
three  amplitude  ratios  for  each  wave  number  ^  '2^*^^  (where  n  has  3  values 
corresponding  to  each  solution  of  equation  A33). 
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4  Write  the  stresses  in  terms  of  particle  displacements,  using  the  constitutive 
relations  and  then  impose  the  stress  free  boundary  condition;  this  results  in  a 
set  of  6  homogeneous  equations. 

5  A  unique  solution  is  available  if  the  determinant  of  the  coefficients  is  equal  to 
zero. 


This  result  along  with  the  equations  of  equilibrium  provides  the  dispersion  rela¬ 
tions  for  the  Lamb  wave  propagation  for  a  fully  anisotropic  plate.  This  case  has  not 
been  fully  solved  due  to  the  complicated  nature  of  the  equations  and  the  lack  of  a 
practical  example.  Special  cases,  such  as  the  case  of  the  monoclinic  crystal,  have  been 
solved  and  detailed  descriptions  are  available. 

The  primary  difference  that  may  be  expected  is  the  existence  of  extra  shear 
modes.  Also  as  has  been  suggested  by  Mindlin  and  Kaul,^^  the  frequency  of  one 
simple  thickness  mode  may  coincide  with  that  of  another  leading  to  difficulties  in 
deriving  the  dispersion  curves.  Nevertheless  work  is  continuing  on  the  orthotropic 
symmetry  case  and  in  particular  the  details  of  the  transverse  resonance  induced  wave 
propagation. 
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